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Abstract 

Starvation is a common stress experienced by bacteria living in natural environments and the ability to adapt to and survive 
intense stress is of paramount importance for any bacterial population. A series of starvation experiments were conducted 
using V. vulnificus 93U204 in phosphate-buffered saline and seawater. The starved population entered the death phase 
during the first week and approximately 1% of cells survived. After that the population entered a long-term stationary 
phase, and could survive for years. Starvation-induced diversification (SID) of phenotypes was observed in starved 
populations and phenotypic variants (PVs) appeared in less than 8 days. The cell density, rather than the population size, 
had a major effect on the extent of SID. SID was also observed in strain YJ016, where it evolved at a faster pace. PVs 
appeared to emerge in a fixed order: PV with reduced motility, PV with reduced proteolytic activity, and PV with reduced 
hemolytic activity. All of the tested PVs had growth advantages in the stationary phase phenotypes and increased fitness 
compared with 93U204 cells in co-culture competition experiments, which indicates that they had adapted to starvation. 
We also found that SID occurred in natural seawater with a salinity of 1%-3%, so this mechanism may facilitate bacterial 
adaptation in natural environments. 
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Introduction 

Vibrio vulnificus is a notorious human pathogen, which causes 
clinical manifestations that range from wound infection to primary 
sepsis, or even death [1]. This species is also a fish pathogen that 
affects eel and tilapia aquaculture in many Asian and European 
countries, including Japan [2], Bangladesh [3], Denmark [4], and 
Spain [5]. It has also been isolated from oysters and shrimps [6,7]. 

However, V. vulnificus is usually found in aquatic environments 
as free-living bacteria. Environmental V. vulnificus isolates are 
reported to be highly heterogeneous. Arias et al. reported high 
genetic diversity after ribotyping 132 strains of V. vulnificus by 
random amplified polymorphic DNA PCR [8]. Wong et al. also 
demonstrated high genetic heterogeneity in environmental isolates 
from the United States and Taiwan using pulsed-field gel 
electrophoresis [9]. The mechanism that allows this high genetic 
diversity to be maintained in marine environments still remains 
unknown. This high genetic heterogeneity may be derived from a 
continuous and abundant supply of mutations in natural 
populations, or it could be the result of relaxed selection in a 
heterogeneous environment. 

In natural environments, bacteria are expected to experience 
repeated "feast or famine" fluctuations, so the capacities to survive 
starvation and respond rapidly to transient nutrient bursts are of 
paramount importance. Gram-negative bacteria react to various 
types of environmental stress with a general stress response. This 
response depends on alternative sigma factors such as RpoS, 
which regulate a suite of genes and redirect resources from 
supporting growth to self-maintenance (reviewed in [10]). These 



include genes related to metabolic adaptations to available 
resources and alterations in the control of motility, cell morphol- 
ogy, stress resistance, and biofilm formation. 

Bacteria activate the general stress response early during 
starvation. They will induce further stringent responses if the 
nutritional limitation persists and the cells will stop to grow, while 
there is a switch to a more error-prone DNA polymerase, which 
may elevate the mutation rate [1 1] . Stress-directed mutation is 
considered a universal mechanism of microbial adaptation to 
environmental stress at the population level during long-term 
starvation [12]. This mechanism may contribute to adaptive 
evolution in Escherichia coli in laboratory culture [13]. This 
phenomenon has been described in many microbes in laboratory 
settings, but its occurrence in natural habitats has not been 
examined previously. 

The population dynamics of V. vulnificus in conditions of 
starvation have been examined in several studies. Biotype 2 strain 
E22 V. vulnificus was reported to survive starvation for 50 days in 
artificial seawater [14]. In a subsequent study, two- to three-order 
drops in population size were seen in E22 after 160-day starvation 
in water with 0.5% or 1.5% salinity at 25°C [15]. Similarly, 
approximately 1% of biotype 1 strain C7184 cells survived at the 
end of 35-day starvation in artificial seawater at 22°C [16]. The 
question is how does starvation affect the physiology of Vibrio spp. 
and does it affect bacterial diversification? Many Vibrio spp. can 
survive long-term starvation, even for several years [15,17]. They 
enter a viable but non-culturable (VBNC) state when they are 
starved at a low temperature [18]. The starved cells are known to 
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change their morphology from rod to coccoid cells, as reported in 
V. vulnificus [15], V. parahaemolyticus [19], V. angustum [20], V. shilloi, 
and V. tasmaniensis [21]. Starved V. parahaemolyticus cells also exhibit 
greater cell adherence and hydrophobicity [9]. These results 
indicate that starvation can lead to dramatic changes in the 
morphology and physiology of Vibrio spp. 

Previously, we isolated V. vulnificus from diseased tilapia [22]. 
The epidemics occurred only in fish maintained in a low salinity 
(<0.5%) environment, which suggests that these tilapia-pathogenic 
isolates were descendants of low salinity-adapted V. vulnificus. 
Further analysis indicated that these isolates formed several genetic 
groups (unpublished data), which indicates either that the ability to 
infect tilapia has evolved more than once, or that horizontal gene 
transfer is responsible for the emergence of novel pathogenic 
strains. We hypothesize that if starvation can lead to phenotypic 
diversification in a natural environment, it may facilitate the 
generation of virulent bacterial variants by providing a larger pool 
of diverse progeny for selection. The aim of the present study was 
to show that starvation can induce the diversification of V. vulnificus 
in laboratory and natural settings, and that this starvation-induced 
diversification (SID) could promote bacterial adaptation and 
persistence. We showed that V. vulnificus survived for an extended 
period and exhibited phenotypic diversification. SID was cell 
density-dependent and a result of adaptive mutation. We provide 
evidence to demonstrate that SID occurred with various salinity 
levels and was observed in two different isolates. 

Materials and Methods 

Bacterial strains, culture media, and growth conditions 

All of the bacterial culture experiments were performed at 
30°C. Bacteria were grown in tryptic soy broth or agar 
supplemented with 1.5% NaCl (TSBS and TSAS, respectively). 
Strain 93U204 is a fish pathogen, which was isolated from 
moribund tilapia collected from Kaohsiung in Taiwan. This isolate 
was a kind gift from Dr. Chia-Ben Chao of Kaohsiung Institute of 
Livestock Disease Control and Prevention, and was identified as V. 
vulnificus by 16S rDNA sequencing (manuscript in preparation). 
YJ016 was a human clinical strain isolated from Tainan, Taiwan, 
and its genome has been sequenced [23]. This strain was kind gift 
from Dr Lien-I Hor of Department of Microbiology and 
Immunology, National Cheng-Kung University, Taiwan. All 
bacterial isolates were properly stored by putting 100 uL overnight 
culture in 1 mL of TSBS supplemented with 50% glycerol, and 
then kept in -80°C freezer (GV039P/M, Kaltis International, 
Taipei, Taiwan). 

Long-term starvation of bacterial cultures 

Bacteria were grown in TSBS with constant 200 rpm horizontal 
shaking (E450, Deng Yng, New Taipei City, Taiwan). Cells were 
harvested in the late exponential phase (approximately 1 x 1 0 9 
CFU/mL), washed twice with phosphate-buffered saline (PBS, 
0.1 M phosphate buffer and 0.85% NaCl, pH 7.3), and re- 
suspended in suitable media. The bacterial cells were diluted in 
PBS to the desired concentrations if needed. 

Most of the starvation experiments were conducted in 30-mL 
Wheaton glass serum bottles (Z 113980, Wheaton, Millville, New 
Jersey, USA). Experiments conducted in other containers are 
indicated. The cells were washed twice with PBS, and re- 
suspended in 10 mL test media. Each serum bottle was capped 
with a flange red rubber stopper (224100-172, Wheaton) and 
sealed with an open aluminum cap (224178-01, Wheaton). The 
cells were incubated at 30°C without shaking and were sampled 
repeatedly using 1-mL insulin syringes (27-gauge, Terumo, Tokyo, 



Japan). The viable counts were determined using the drop plate 
method [24]. The bacterial suspensions were serially diluted in 
TSBS and five 10 U.L drops of each dilution were placed on TSAS 
plates. The plates were incubated at 30°C for 18 h before they 
were counted. We chose to use TSBS as diluent to prevent possible 
loss of viability of these stressed cells on agar plate, which has been 
documented for V. vulnificus [25]. The serial dilution process 
generally took less than 10 minutes to perform, and in a 
preliminary test, cells starved in PBS did not increase significantly 
(n = 6, paired t test, P>0.05) after being transferred into TSBS for 
20 minutes. Therefore we consider using TSBS as diluent did not 
invalidate the drop plate results. To compare the bacterial survival 
in various media, 93U204 cells were re-suspended in PBS, natural 
seawater with 3% salinity (NSW), half-strength natural seawater 
diluted with distilled water (HNSW), or TSBS. The salinity of 
seawater was determined using a hand-held salinity refractometer 
(MR100ATC, Milwaukee Instruments, Rocky Mount, North 
Carolina, USA). The starting concentration of 93U204 was 
approximately 5xl0 9 CFU/mL. The survival rate was deter- 
mined every day for 7 days by plate counts on TSAS after serial 
dilution. In the long-term starvation experiment, we re-suspended 
the cells in PBS or TSBS for 214 days using a starting 
concentration of lxlO 9 CFU/mL. The bacterial viability was 
measured on days 0, 1, 2, 3, 8, 14, 35, 56, 93, 115, and 214. 

Phenotypic characterization 

We developed a Motility-Hemolysis-Proteolysis (MHP) group- 
ing system to classify V. vulnificus isolates, which was based on three 
readily identifiable characteristics, namely motility (M), hemolysis 
(H), and proteolysis (P). The isolates were classified into one of 100 
phenotype groups according to their motility (five levels, from M0 
to M4), hemolysis (four levels, from HO to H3), and proteolysis 
(five levels, from P0 to P4) activity levels. Motility, hemolytic 
activity and proteolytic activity were measured on motility agar 
(1% tryptone, 2% NaCl, 0.25% agar), blood agar (LB agar 
supplemented with 5% sheep blood), and skim milk agar (LB agar 
supplemented with 2% skim milk). In every test the progenitor 
93U204 was inoculated on the same test plate as a control, and the 
diameters of its colony on motility agar, clear zone on blood agar 
and clear zone on skim milk agar were recorded. We defined M0 
to M4 as having colony size of 0, 0-0.33, 0.33-0.67, 0.67-1.2, or 
over 1.2 times of that measured in 93U204 control; HO to H3 as 
having clear zone on blood agar of 0, 0-0.5, 0.5-1.5, or over 1.5 
times of that measured in 93U204 control; and P0 to P4 as having 
clear zone on skim milk agar of 0, 0-0.33, 0.33-0.67, 0.67—1.2, or 
over 1.2 times of that measured in 93U204 control. The 
progenitor 93U204 had the M3H2P3 phenotype, and the human 
pathogenic strain YJ016 had the M3H1P3 phenotype. Some 
93U204 phenotypic variants (PVs) developed opaque or translu- 
cent colonies. These PVs included an additional "O" for opaque 
and "T" for transparent in their MHP label. 

To determine the phenotypic composition, colonies from each 
sample were picked randomly with sterile toothpicks and 
inoculated into 150 uL LB in a 96-well microtiter plate, followed 
by overnight incubation at 30°C. We used a stainless steel 96-pin 
replicator to transfer the bacterial cells onto motility agar, skim 
milk agar, blood agar, and thiosulfate-citrate-bile salts-sucrose 
(TCBS) agar. The results for motility, proteolytic activity, 
hemolytic activity, and growth on TCBS were determined at 6, 
20, 48, and 20 h after inoculation at 30°C. Isolates with the same 
MHP grouping were defined as a PV. 
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Starvation-induced diversification (SID) experiments 

In this study six different SID experiments in PBS were 
conducted. Each bottle of inoculated PBS was treated as a 
population. Populations within the same experimental group 
received the same 93U204 preparation as the inocula. For each 
population, a sampling series was conducted at the indicated time 
for each experiment to trace the dynamics of phenotypic change. 
For each sample, approximately 30 colonies were isolated for 
phenotypic characterization and the result was used to calculate 
the relative abundance of each PV. The starting bacterial density 
was adjusted to approximately 10 9 CFU/mL, and the experiments 
were conducted in 10 mL of PBS in 30-mL serum botdes, unless 
indicated otherwise. 

The settings for each experiment were as follows. (1) The 
analysis of temporal change in the PV composition included two 
sets of experiments, with one population in each set. Both 
populations received the same inoculum. Samples were taken on 
days 0, 8, 16, 24, and 32. (2) The analysis of SID in replicate 
populations comprised two sets of experiments, with five 
populations in each set. Samples were taken on day 33. (3) The 
SID and population survival rate analysis comprised 12 sets of 
experiments, with 1-3 populations in each set and a total of 20. 
Samples were taken on days 14 or 15. (4) The analysis of the effect 
of cell density on SID comprised four sets of experiments, with 
three populations in each set. The bacterial densities in these four 
sets were adjusted to approximately 10 9 , 10 7 , 10', and 10 3 CFU/ 
mL. Samples were taken on days 14, 28, and 42. (5) The analysis 
of the effect of population size on SID comprised three sets of 
experiments, with three populations in each set. The bacterial 
suspension volumes used in the three sets were 100 mL, 10 mL, 
and 3 mL, in 125-mL serum botdes (Zl 14014, Wheaton), 30-mL 
serum bottles, and 30-mL serum bottles, respectively. Samples 
were taken on days 14, 28 and 42. (6) Two V. vulnificus strains, 
93U204 and YJ016, were used to test the effect of salinity on SID. 
For each strain, the analysis comprised one set of experiments, 
with three treatments in each set. Samples were taken on days 14, 
28, and 42. These sets used natural seawater with 1%, 2%, and 
3% salinities. The seawater was sterilized using 0.22 urn filter. The 
salinity was checked using a refractometer. The desired salinity 
was achieved by mixing freshly collected seawater with freshwater 
from the same estuary area. 

Competition experiments 

Three isolates, each from a separate PBS-starvation survivor 
population of 93U204, were used in the competition experiments. 
B33-12 was isolated from the trial B 33-day population; C30-2 was 
isolated from the trial C 30-day population; and D33-4 was 
isolated from the trial D 33-day population. All of the strains 
derived from the starvation cultures were identified as V. vulnificus 
by species-specific vvp PCR [26]. Another primer pair of pi 34/ 
pi 35 was also used, which targeted V. vulnificus rpoS [27]. The rpoS 
primer pair comprised the forward primer pl34 (5'-ACATAAG- 
gATAATTACCTCAgTgC-3 '), which started 264 bp upstream 
of the YJ016 rpoS, and the reverse primer pl35 (5'- 
CgTCAAAAATTgTCTCAACTTg-3'), which ended 276 bp 
downstream of YJ0 1 6 rpoS and was located in mutS. The expected 
size of this PCR product was 1530 bp. The PCR protocol 
comprised a preheating step at 94°C for 6 min, followed by 30 
cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 2 min, with a 
final extension at 72°C for 7 min. Their rpoSDNA sequences were 
identical to those of the 93U204 progenitor cell according to 
nucleotide sequencing with an Applied Biosystems 3730x1 DNA 
Analyzer, which was conducted at the VYM Genome Research 
Center of National Yang-Ming University. The obtained nucle- 



otide sequences were compared with the GenBank DNA database 
using the BLASTn program. 

All of the isolates were stored at — 80°C in TSBS with 15% 
glycerol. The stock was subcultured in 12 mL TSBS in a 125-mL 
flask and grown to the late exponential phase with 200 rpm 
horizontal shaking. The culture was then washed twice with PBS 
by centrifugation (10 min, room temperature, 4500 Xg). The 
bacterial preparation was re-suspended in PBS to the desired 
concentration. The progenitor 93U204 was mixed at 1000:1 with 
each of the isogenic cultures derived from C30-2, D33-4, or B33- 
1 2 in PBS at the start of a two-week experiment. Monocultures of 
the 93U204 progenitor and mutants were included as controls. 

The bacterial concentration of each sample was determined 
using the drop plate method with TSAS plates. The relative 
abundance of each competing population was determined every 2 
days for 14 days. We used phenotypic characterization to identify 
the relative amount of the two competing phenotypes, because 
only a small fraction of the cells changed their phenotype up to the 
late stage of the starvation experiment. In each of these 
competition experiments, cells with reduced motility and en- 
hanced proteolytic activity were considered to be descendents of 
C30-2, cells with reduced motility that were unable to grow on 
TCBS were considered to be descendents of D33-4, and cells with 
reduced motility and enhanced proteolytic activities were consid- 
ered to be descendents of B33-12. 

Diversity evaluation and statistical analysis 

The bacterial diversity was evaluated using the Shannon 
diversity index, H' = —^Jj^np,, where pi is the relative abundance 
of each phenotype. Each phenotype was defined as having the 
same motility, hemolysis, and proteolysis activity characteristics. 
Cluster analysis and linear regression were performed using PAST 
software (http://folk.uio.no/ohammer/past/) [28]. 

Electron microscopy 

All of the electron microscopy experiments were conducted in 
the Electron Microscopy Center, Tzu-Chi University. A bacterial 
suspension from each preparation was centrifuged at 5,000 for 5 
min at 4°C to collect bacterial cells. The pellets were re-suspended 
and fixed with 2.5% glutaldehyde in PBS at 4°C. On the day the 
observations were made, the bacterial cells were washed twice with 
0.1 M phosphate buffer with 5% sucrose, post-fixed with 1% 
osmium tetraoxide for 1 h, and washed twice with 0.1 M 
phosphate buffer with 5% sucrose. Next, a 5 |jL drop of bacterial 
suspension was placed on a Formvar/ carbon-coated grid for 2 min 
to allow absorption of bacteria onto the grid. A drop of 2% 
phosphotungstic acid was placed on the grid to achieve negative 
staining. The grids were then air-dried and examined using a 
Hitachi H-7500 transmission electron microscope (Hitachi, 
Tokyo, Japan) at 80 kV. Digital images were collected using a 
2048x2048 Macrofire monochrome CCD camera (Optronics, 
Goleta, California, USA). Cell diameter and length measurements 
were made using these images with ImageJ v. 1.46 (http:/ /rsbweb. 
nih.gov/ ij/index.html). 

Chemicals and reagents 

All of the PCR-related reagents were included in the Taq DNA 
Polymerase 2X Master Mix Red kit (1.5 mM MgCl 2 , Ampliqon 
IIII, Odense, Denmark). The electron microscopy chemicals and 
reagents were obtained from Electron Microscopy Services 
(Hatfield, Pennsylvania, USA). The culture media and ingredients 
were obtained from Becton Dickinson and Company (Sparks, 
Maryland, USA), while all other chemicals were obtained from 
Sigma (St Louis, Missouri, USA). All of the primers were 
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purchased from Gene Messenger (Kaohsiung City, Taiwan). Fresh 
sheep blood was obtained from Creative Media Products Ltd (New 
Taipei City, Taiwan). 

Results 

Survival of V. vulnificus 93U204 after starvation in various 
culture media 

The repeated appearance of new V. vulnificus pathogens in 
tilapia farms and their high heterogeneity in natural populations 
(unpublished data) suggests that some environmental biotype 1 
strains might repeatedly invade and persist in low-salinity tilapia 
farms, which later became pathogens. Thus, we aimed to 
determine whether environmental stresses, such as starvation, 
could contribute to the emergence of virulence. 

First, we examined the survival of V. vulnificus 93U204 in 
nutrient-poor PBS and seawater, or nutrient-rich TSBS. The 
number of culturable 93U204 in PBS dropped rapidly by two 
orders from the starting concentration of 5 x 1 0 9 CFU/ mL during 
the first week in PBS, NSW and HNSW (Figure 1A). The viable 
population was even smaller in TSBS, to approximately 8xl0 2 
CFU/mL on day 1 and below the detection limit (20 CFU/mL) 
after day 3. 

To understand the long-term survival of V. vulnificus, we 
inoculated PBS and TSBS with cells at approximately lxlO 9 
CFU/mL, and tracked the viable population dynamics for 7 
months (Figure IB). Three-order and five-order drops were seen 
after 7 days in the PBS and TSBS groups, respectively. The 
mortality stabilized after 2 weeks, then the populations entered a 
slow and gradual declining phase in both groups. In another 
experiment with a starting cell density of 5 x 1 0 8 CFU/ mL, 
approximately 10 CFU/mL of V. vulnificus survived starvation in 
PBS at 1215 days after inoculation. We conclude that a starved 
93U204 population will typically experience a rapid decrease 
during the first week, before the mortality stabilizes after the 
second week. The population then enters a slow period of decline, 
but can survive for months to years. 

We examined the morphology of starved cells by transmission 
electron microscopy (Figure 1C). The un-stressed log phase cells 
exhibited the typical vibrioid morphology, whereas the cells 
starved in PBS for 4 days became more slender. Most of the cells 
from the 4 day TSBS population were ruptured, which indicated 
their extremely low viability in this treatment. PBS-starved cells 
were a similar length to the log-phase cells, but their cell diameters 
were much smaller (Figure ID). 

Starvation and PBS-induced phenotypic diversification in 
V. vulnificus 

The high mortality in starved V. vulnificus populations may exert 
strong selection on survivors. We conducted 2 independent PBS 
starvation trials of 93U204 and monitored the population for 
phenotypic change. SID was seen in both trials (Figure 2A). A total 
of 16 and 13 PVs were found in trial 1 and trial 2, respectively. We 
designated the progenitor 93U204 M3H2P3 phenotype as the 
"wild-type (WT)" phenotype. PV was first detected at day 8 in trial 
1, and at day 16 in trial 2. The proportion of WT cells decreased 
with time as the number of detected PVs increased. To confirm 
that these changes were not transient, we picked five morpholog- 
ically distinct colonies from each population to assess their 
phenotypic consistency through subcultures. After at least three 
subcultures onto TSAS, none reverted back to the WT phenotype, 
which demonstrated that these changes were not transient and 
were likely to be genetic. 



The PVs in each population were categorized further to detect 
general trends in the motility, hemolysis, or proteolysis activity 
(Figure 2B). PVs with reduced motility (mot-R PVs) and with 
reduced hemolytic activity (hem-R PVs) increased with time, and 
the mot-R PVs appeared earlier than the hem-R PVs during 
starvation. PVs with enhanced proteolysis activity (pro-E PVs) 
appeared transiently on day 24 and disappeared or decreased 
greatly on day 32 in both trials. We conclude that a reduction in 
motility, a reduction in hemolytic activity, and an enhancement in 
proteolytic activity could promote the survival of 93U204 during 
starvation. 

SID was associated with cell survival in populations 

If diversification is a bacterial strategy for overcoming 
environmental stress, there should be a correlation between 
diversity and population survival. Thus, we collected survival and 
PV composition data from 20 starved 93U204 populations. All of 
these populations were started with approximately 10 9 CFU/mL 
in 10 mL PBS, followed by starvation for 14 or 15 days before 
sampling. Three of the 20 populations had very low viability at 
approximately 0.1% of the original population (Figure 2C). SID 
was seen in all but these three populations, which remained 
homogeneous and all of the colonies examined (30 for each 
population) had an unchanged phenotype. This result suggests that 
SID may be needed for the population to overcome environmental 
stress. 

Source of PVs in starvation survival populations 

PVs found in starved population may have been present in low 
numbers in the original population, or they may have been 
mutants derived from WT cells during starvation. Thus, we 
conducted five parallel starvation experiments using the same 
original population as the inoculum (Figure 3). We hypothesized 
that if these PVs were original members of the inoculum 
population and were selected during starvation, the compositions 
of these five surviving populations should be very similar. 

The results showed that each population had a unique 
phenotypic composition after 33 days of starvation (Figure 3B), 
which suggested that these PVs were novel mutants that evolved 
independendy in each population. However, similar trends were 
seen among the populations (Figure 3C). In both trials, mot-R PVs 
dominated the five populations, which indicates that a reduction in 
motility may provide more advantages than other mutations 
during early starvation. In some populations, hem-R PVs also 
increased in proportion (Figure 3B). 

Advantages of starvation survivors compared with wild- 
type cells 

Certain PVs were observed repeatedly in the starvation 
experiments (Figures 2 and 3), which suggests that these PVs 
may have been more resistant to starvation-induced death, or that 
they may have been related to traits that made them more 
competitive than WT cells during starvation. To demonstrate that 
PVs performed better than the WT strain during starvation 
conditions, we selected three isolates, i.e., C30-2, D33-4, and B33- 
12, from three different survivor populations and performed 
competition tests with WT cells. C30-2 was isolated from the trial 
C 30-day population, and had reduced motility and enhanced 
proteolytic activity (M1H2P4 phenotype). D33-4 was isolated from 
the trial D 33-day population, and had reduced motility, reduced 
hemolysis, and increased proteolytic activity, and could not grow 
on TCBS plates (M1H1P4 phenotype). B33-12 was isolated from 
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Figure 1. Survival of V. vulnificus cultures in various media. (A) Dynamics of V. vulnificus populations in phosphate-buffered saline (PBS, •), 
natural seawater (NSW, <>), half-strength natural seawater (HNSW, □), and tryptic soy broth with 1.5% NaCI (TSBS, A) during the first week after 
inoculation. The results shown were derived from three replicate populations. (B) Dynamics of V. vulnificus populations in PBS (•) and TSBS (A) during 
an extended period of starvation (214 days). Representative results are shown from one experiment. (C) Electron micrographs of V. vulnificus cells 
showing the cellular morphology of a log-phase population (Log phase), a population starved in PBS for 4 days (PBS), and a population starved in 
TSBS for 4 days (TSBS). Most of the cells starved in TSBS were ruptured. All of the bars represent 2 um. (D) The diameters and lengths of log-phase 
cells (A) and cells starved in PBS for 4 days (•). 
doi:1 0.1 371 /journal.pone.0088658.g001 



the trial B 33-day population, and had reduced motility and 
hemolysis activity levels (M1H1P3 phenotype). 

When cultured alone in PBS, the survival rates of all the tested 
PVs were similar to that of the progenitor control (Figure 4, 
upper), which indicated that they were not more resistant to 
starvation-induced death. When inoculated at 0.1 % in a mixed 
culture with progenitor 93U204, the C30-2-like M1H2P4 cells 
exhibited a slow increase and reached 14% on day 14 (Figure 4A) 
(95% and 81% in the second and third trials), while D33-4-like 
M1H2P4 cells reached 56% on day 14 in mixed culture (Figure 
4B) (79% and 62% in the second and third trials, respectively). 
B33-12-like M1H1P3 cells reached 40% on day 2 and continued 
to increase until they reached 95% on day 14 (Figure 4C) (49% 
and 79% in the second and third trials, respectively). The relative 
abundance of each mutant varied among trials, but they 
dominated from 0. 1 % in all three trails. In mutant monoculture, 



no isolates switched back to the WT phenotype. Similarly, only 
less than 6.5% of mutant-like isolates were detected in the 93U204 
progenitor monoculture at the end of the experiments. Therefore, 
we conclude that all three tested PVs were more competitive than 
the WT cells at surviving starvation. 

Cell density, but not the number of cells, had a major 
effect on SID 

High density can be stressful for bacterial cells and lead to more 
mutations in the population. Thus, we tested the effect of cell 
density on the generation of SID. We examined SID at 10 , 10 7 , 
10 5 , and 10 3 CFU/mL. The PV diversity appeared to increase 
with density (Figure 5A). The Shannon index correlated positively 
with log-transformed cell density (linear regression, r= 0.76684, 
P = 0.0048). Both percentage of WT cells (r =-0.81659, 
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Figure 2. Phenotypic diversification of the V. vulnificus 93U204 
population during starvation in phosphate-buffered saline. (A) 

Relative abundance of phenotypic variants in starved populations from 
two replicate trials. The population gradually diversified from the 
93U204 M3H2P3 phenotype. (B) Compositions of populations, which 
were categorized based on their motility (top panel), hemolytic activity 
(middle), and proteolytic activity (bottom panel). WT, phenotype of the 
parental 93U204. (C) Association of diversity with mortality in starved 
populations. Twenty populations were starved in phosphate-buffered 
saline, and their survival and diversity were measured on days 14 or 15. 
The figure shows the distinction between low survival/un-diversified 
populations (left) and high survival/diversified populations. These 20 
populations were derived from 12 separate inoculations (indicated by 
different colors). 
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P = 0.0008) and ratio of final to initial population size (r = — 
0.88628, P = 0.0001) correlated negatively with log-transformed 
cell density. However, this trend was only evident in the day 14 
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Figure 3. Phenotypic variant (PV) compositions of different 
sets of starvation populations. (A) The PV compositions of two sets 
of starvation populations, which started with the same inoculum, were 
analyzed for 33 days (n = 5). (B) Summaries of the phenotypic 
composition in terms of the motility (M), hemolysis (H), and proteolysis 
(P) activity levels are shown in. (C) Cluster analysis of the PV 
composition of each population. Green rectagulars, populations 
received inoculum 1; organge, populations received inoculum 2; blue, 
unstarved control. The analysis was based on the Bray-Curtis similarity 
and paired group methods. 
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population and was absent on days 28 and 42, which suggests that 
the density effects were probably masked by saturation because of 
continuous diversification. 

In these experiments, the increase in density was confounded by 
the increased absolute number of cells in the populations. Larger 
populations contained more cells so they were more likely to 
contain more mutants, which meant that there was a higher 
likelihood of beneficial mutations accumulating in these popula- 
tions. In the next experiment, we fixed the starting cell density at 
10 9 CFU/mL and varied the volume of the medium among 100, 
1 0, or 3 mL, before examining the effect of the population size on 
SID (Figure 5B). Neither Shannon index nor percentage of WT 
cell correlated significantly with medium volume. Therefore, we 
conclude that the cell density has a more important effect than the 
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population size during phenotypic diversification in conditions of 
starvation. 

SID in seawater 

SID was readily observed in a laboratory setting, so we 
investigated whether this phenomenon also occurs in natural 
seawater. We starved V. vulnificus 93U204 cells in natural seawater 
samples, which were adjusted to salinities of 1%, 2%, or 3%, and 
we monitored any phenotypic changes (Figure 6A and B). The 
SID observed in natural seawater was comparable to that in PBS, 
which indicates that SID can occur in natural marine environ- 
ments. For comparison, another V. vulnificus isolate YJ016 was also 
used in this starvation experiment (Figure 6C and D). SID was 
observed in both 93U204 and YJ016, although YJ016 diversified 
at a faster rate. The three populations cultured with different 
salinities had similar SID, which suggests that this phenomenon 
was not affected by salinity. 

Discussion 

In this study, we determined the population dynamics of 
93U204 during long-term starvation. The survival rates were 
similar to those reported in previous studies [14—16]. During the 
first week, there was a dramatic decrease in the viable population. 
Two weeks after inoculation, the population increased slightly, 
before the population entered a long and gradual decline (Figure 
1). The life cycle of E. coli has been classified into lag, log, 
stationary, death, and long-term stationary phases [29] . E. coli cells 
grow rapidly during the log phase but the growth slows in the 
stationary phase when the nutrient supply decreases, and finally 
the cells progress to the death phase. Typically, >99% of the 
population lose their viability during the death phase, but the 
remaining cells can survive for several years with few declines in 
their cell numbers. The survival of V. vulnificus followed similar 
dynamics to those previously reported for E. coli, with an apparent 
death phase (rapid decrease) and a long-term stationary phase 
(slow decrease). A starved E. coli population maintained its size at 
approximately 1 0 6 CFU/ mL for over 5 years when sterile distilled 
water was supplied to maintain the volume and osmolarity [29]. In 
our experiments, V. vulnificus remained at approximately 10 6 
CFU/mL in PBS and 10 5 CFU/mL in TSBS for 214 days (Figure 
IB), while it remained at approximately 10 CFU/mL in PBS for 
1215 days. Therefore both species have comparable survival in 
longterm starvation. 

After treating the stationary phase E. coli population with 
antibiotics to remove any dividing cells, some non-dividing cells 
were found to be dividing [30]. These cells were named 
"persisters" that allowed the population to survive in a harsh 
environment. Similar persisters may also play a role in the earlier 
stages during the survival of V. vulnificus, but they are not likely to 
be the major source of survivors because the "revived" persisters 
had the WT phenotype. The divergence between parallel 
populations and the fixed phenotype during passages strongly 
suggests that PVs are mutants derived from the progenitor cells. 
The cells acquired adaptive mutations during starvation, before 
subsequent selection, and they survived the stressful conditions. 

Vibrios are known to enter the VBNC state under stress. In V. 
vulnificus the VBNC state could be induced through a temperature 
downshift to 5°C [31] or 4°C [32]. Although the cell viability seen 
in our starvation experiment may be similar to that in VBNC 
population, it is unlikely that the loss of viability is due to VBNC. 
Our experiments were conducted at 30°C, and no VBNC V. 
vulnificus has been reported at this temperature so far. An elevation 
of temperature to 21°C is sufficient to resuscitate VBNC V. 



vulnificus within 24 h [33]. We compared the property of 93U204 
cells kept in PBS at either 30°C or at 4°C. After 9-day starvation 
treatment, most cells kept at 30°C co-stained with both propidium 
iodide and SYTO 9 (using LIVE/DEAD fecLight Bacterial 
Viability Kit, for microscopy & quantitative assays, L7012, Life 
Technologies, Carlsbad, California, USA)(unpublished data), but 
co-staining was rarely seen in cells kept at 4°C. Addition of 
catalase is known to increase the survival of VBNC cells on agar 
plates [34] . Addition of catalase increased survival of cells kept at 
4°C (with vs. without catalase, paired t-test, P = 0.038) by 4.6 ±3.2 
times, but it had no effect on cells kept at 30°C (paired t-test, 
P = 0.4655). These data suggest that the loss of viability seen in our 
starvation experiments was not due to the induction of VBNC 
state. 

Most cells that survived the death phase exhibited phenotypic 
changes (Figure 2). These changes were generally advantageous 
compared with the progenitor cells and they are referred to 
collectively as growth advantage in the stationary phase (GASP) 
phenotypes [29]. Cells that exhibit GASP phenotypes during 
starvation have been reported in bacteria such as Campylobacter 
jejuni [35], Pseudomonas [36,37], Geobacter [38], and enterobacteria 
[39]. A GASP phenotype was also observed in V.fischeri [40]. A 
non-toxigenic V. cholerae isolate derived from a starved long-term 
slab culture was shown to have growth advantages compared with 
its toxigenic progenitors [41]. Mutations that lead to GASP have 
been studied in E. coli (reviewed in [29]). The loss of a functional 
rpoS may be advantageous in increasing amino acid utilization, 
while hp and the ybeJ-gltJKL cluster encode the leucine-responsive 
protein and a high-affinity aspartate and glutamate transporter, 
respectively. In Vibrio spp., however, there is still no known genetic 
basis for GASP phenotypes. 

Three types of strategies may facilitate bacterial survival during 
starvation: persister, scavenger, and replicator strategies. Cells that 
employ a persister strategy can persist during starvation conditions 
via physiological adaptation. These cells may have to trade their 
nutrient acquisition ability (nutritional competence) for self- 
preservation, according to the SPANC balance hypothesis [42]. 
Scavengers are cells that can increase their nutrient input, either 
by acquiring the ability to explore novel resources or by increasing 
their substrate-binding affinity or substrate-transporting ability. 
Replicators are cells that reduce their energy expenditure via 
mutation in functions that are not required for survival so they 
may have higher rates of replication. Scavengers and replicators 
may replicate faster than their progenitor cells, so they are more 
likely to dominate the surviving population. 

Several systems have been used to study bacterial responses to 
nutrient limitation. Lenski's laboratory have conducted a long- 
term evolution experiment (LTEE) with E. coli [43] . In this system, 
the bacteria are transferred to new broth each day. During the first 
few hours, replicators may have advantages, but the scavengers are 
favored over time. Persisters have an advantage only when the 
population experiences severe nutrient deprivation, which is not 
likely to happen in the LTEE. Therefore, replicators and 
scavengers are expected to dominate in this system. The mean 
fitness increased by approximately 35% in the 2000-generation 
population [44] . The adaptation led to a decrease in the lag time 
before growth and an increase in the maximum growth rate, 
although the death rate and the nutrient concentration required to 
support half the maximum rate did not change significantly. 
Therefore, replicators were favored in this system. 

Ferenci's laboratory pioneered the use of a chemostat system to 
study the response and evolution of E. coli in a constant low- 
nutrient environment [45]. The nutrient concentrations were 
maintained at stable low levels, so the scavengers had an 
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Figure 6. Starvation-induced diversification (SID) of V. vulnificus strains 93U204 and YJ016 in seawater. We tested the SID in natural 
seawater samples with salinities of 1%, 2%, or 3%. The results for 93U204 are shown in (A) and (B), while those for YJ016 are shown in (C) and (D). The 
starting concentration was 10 9 CFU/mL. The phenotypic variant compositions (A and C) and the summarized population compositions (B and D) for 
motility (upper), hemolysis (middle), and proteolysis activity (lower) are also shown. 
doi:1 0.1 371 /journal.pone.0088658.g006 



advantage over replicators and persisters. In this system, an rpoS 
mutant rapidly became dominant within the first 3 days [46]. 
Growth in these conditions is highly r/xiD-dependent, which 
supports Commutations [47]. Mutations in the rpoS, mlc, malT, and 
mgl have been demonstrated to increase the expression of ompF, 
lamB, mglBAC, and ptsG, which improves high-affinity glucose 
uptake [48 _ 50]. Thus, the scavengers dominated this system. 

In our study, we starved the cells in PBS from the start, so 
persisters may have possessed advantages compared with other 
cells during the first few days. Dead cells accumulated rapidly over 
time and the surviving cells actively utilized the nutrients released 
from the dead cells, so they gradually increased in number and 
dominated the population. The first types of PVs that appeared in 
the starved populations were mot-R PVs (Figures 2 and 6). 
Reductions or loss of motility after starvation have been reported 
in diverse bacteria, including Rhizobium meliloti [51], V. angustum 
[52], and unidentified marine bacteria [53]. Flagellar motility 
requires much energy, so a reduction in bacterial motility has the 
benefit of saving energy, which is consistent with the recendy 
proposed Black Queen hypotheses [54]. However, these PVs did 
not have an increased growth rate in TSBS. It is likely that these 
PVs lacked a growth advantage, or they may only possess a growth 
advantage in low-nutrient conditions. During a stringent response, 
E. coli downregulates the expression of flagellar genes [55], [56]. 
Recently, it was reported that V. vulnificus downregulates its 
flagellar regulator FlhF in the stationary phase via the quorum 
sensing master regulator SmcR [57]. This indicates that cells are 
not motile in the stationary phase so those that lose motility will 
not be selected against and will tend to be detected. Without 
appropriate evidence, we were not able to tell whether the rapid 
appearance of mot-R PVs was a result of a replicator or scavenger 
strategy. However, based on the presence of large numbers of dead 
cells, a scavenger strategy was more favored. 

In the present study, the average proteolytic activity of the 
93U204 survivor populations remained unchanged or increased, 
but the hemolytic activity decreased with time (Figure 2). These 
trends may represent a trade-off between using energy to produce 
the corresponding proteins and acquiring nutrients from dead cells 
with these proteins. In the set of experiment shown in Figure 6B, 
the proteolytic activity and hemolytic activity of the survivors 
remained largely unchanged, although mot-R PVs still dominated. 
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